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^^ , Abstract. The Spitzer Space Telescope is pushing into new frontiers in high 

^^ ■ redshift astronomy. Closer to home, Spitzer is making an equally large impact 

Cn ' on our understanding of galaxy formation and evolution. In this contribution we 

$H , present mid-infrared diagnostics based largely on data from the Spitzer Infrared 

Q-e Nearby Galaxies Survey (SINGS). Our main result is that these mid- infrared 

•^T diagnostics effectively constrain a target's dominant power source. The combi- 

nation of a high ionization line index and PAH strength serves as an efficient 
discriminant between AGN and star-forming nuclei, confirming progress made 
with ISO spectroscopy on starbursting and ultraluminous infrared galaxies. The 
sensitivity of Spitzer allows us to probe fainter nuclei and star-forming regions 
qq • within galaxy disks. We find that both star-forming nuclei and extranuclear 

f-^ | regions stand apart from nuclei that are powered by Seyfert or LINER activ- 

f~^ , ity. In fact, we identify areas within three diagnostic diagrams containing >90% 

t^J- ■ Seyfert/LINER nuclei or >90% H II regions/H II nuclei. We also find that, com- 

f^ ' pared to starbursting nuclei, extranuclear regions typically separate even further 

VO . from AGN, especially for low-metallicity extranuclear environments. In addition, 

^^ ■ instead of the traditional mid-infrared approach to differentiating between AGN 

and star- forming sources that utilizes relatively weak high-ionization lines, we 
show that strong low-ionization cooling lines of X-ray dominated regions like 
[Si II] 34.82 /im can alternatively be used as excellent discrimants. 
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Introduction 



The goal of this study is to explore whether mid-infrared diagnostics developed 
for luminous/ultraluminous infrared galaxies and bright Galactic H II regions 
can be improved upon and extended to the nuclear and extra-nuclear regions 
within normal and infrared-faint galaxies. A traditional method for character- 
izing a galaxy's nuclear power source uses ratios of optical emission lines such 
as [O II] 3727A, H/3 486lA, [O III] 5007A, [O I] 6300A, Ha 6563A, [N II] 6584A 
and [S II] 6716,673lA (e.g., Baldwin, Phillips, & Terlevich 1981; Veilleux & Os- 
terbrock 1987; Ho, Filippenko, & Sargent 1997; Kewley et al. 2001; Kauffmann 
et al. 2003). A plot of [O III] /H/3 versus [N II] /Ha, for example, will typically 
separate Seyferts, LINERs, and starburst nuclei. Since nuclei are often heavily 
enshrouded by dust, especially in luminous and ultraluminous infrared galaxies 
(LIRGs and ULIRGs), an important limitation to galaxy optical diagnostics is 
the effect of extinction. In anticipation of the data stream from space-based in- 
frared platforms, early theoretical work with photoionization models showed that 
infrared ionic fine structure line ratios could profitably enable astronomers to ap- 
proach galaxy classification from a new perspective (e.g., Voit 1992; Spinoglio & 
Malkan 1992). The advent of sensitive infrared line data from the Infrared Space 
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Observatory was an important first step to peering more deeply into buried nu- 
clear sources (Genzel et al. 1998; Laurent et al. 2000; Sturm et al. 2002; 
Peeters, Spoon, & Tielens 2004). Genzel and collaborators were the first to 
show that ionization-sensitive indices based on mid-infrared line ratios correlate 
with the strength of polycyclic aromatic hydrocarbon (PAH) emission features. 
AGN in particular show weak PAH and large ratios of high-to-low ionization line 
emission. High ionization lines are not the only route, though, to determining 
whether a galaxy harbors a strong AGN. Similar to how the [O I] 6300A line is 
an AGN diagnostic, the [O I] 63 /im far-infrared line can also help to decipher 
a galaxy's power source (Dale et al. 2004a). Mid- infrared lines observed by 
ISO were also used to probe the physical characteristics and evolution of purely 
starbursting nuclei (Thornley et al. 2000; Verma et al. 2003) and Galactic H II 
regions (Vermeij &: van der Hulst 2002; Giveon et al. 2002). An important result 
stemming from these efforts is that stellar aging effects appear to result in H II 
regions generally having higher excitations than starbursting nuclei. 

The unprecedented sensitivity and angular resolution afforded by the Spitzer 
Space Telescope allow even more detailed studies of galaxy nuclei (e.g., Armus 
et al. 2004; Smith et al. 2004). SINGS takes full advantage of Spitzer's ca- 
pabilities by executing a comprehensive, multi-wavelength survey of 75 nearby 
galaxies spanning a wide range of morphologies, metallicities, luminosities, and 
star formation activity levels (Kennicutt et al. 2003). The sensitivity of Spitzer 
coupled with the proximity of the SINGS sample allows dwarf galaxy systems 
fainter than Lfir ~ 10 7 L & to be spectroscopically probed in the infrared for the 
first time. In addition, prior to Spitzer the only individual extragalactic H II re- 
gions that were detectable with infrared spectroscopy resided in the Local Group 
(e.g., Giveon et al. 2002; Vermeij et al. 2002). In contrast, SINGS provides in- 
frared spectroscopic data for nearly 100 extragalactic H II regions, residing in 
systems as near as Local Group members to galaxies as far as ~25 Mpc. The 
SINGS dataset thus samples a wider range of environments than previously ob- 
served with infrared spectroscopy. This diversity in the SINGS sample provides 
a huge range for exploring physical parameters with the mid-infrared diagnostics 
presented below. 



2. The Sample and Data 

The sample of nuclear targets analyzed in this study derive from the SINGS 
Third Data Release. These 50 nuclei come from a wide range of environments 
and galaxies: low-metallicity dwarfs; quiescent ellipticals; dusty grand design 
spirals; Seyferts, LINERs, and starbursting nuclei of normal galaxies; and sys- 
tems within the Local and M 81 groups (see Kennicutt et al. 2003). The 26 
extranuclear sources studied in this work also come from the SINGS Third Data 
Release. These targets stem from the original set of 39 optically-selected sources 
listed by Kennicutt et al. (2003). The optically-selected OB/H II regions cover 
a large range of metallicity, extinction-corrected ionizing luminosity, extinction, 
radiation field intensity, ionizing stellar temperature, and local H2/H I ratio as 
inferred from CO. Further details on these data and the archival data used in 
this study are fully described in Dale et al. (2006). 
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2.1. Mid-Infrared Spectral Diagnostics 

Emission Line Ratios and PAH Strength Studies show that PAH features are 
quite prominent throughout much of the interstellar medium except for regions 
characterized by exceptionally hard radiation fields such as those that arise in 
AGN and the cores of H II regions (Cesarsky et al. 1996; Sturm et al. 2000). 
Furthermore, the different energies required to stimulate [O IV]25.89 /jm and 
[Ne II] 12.81 /im means their flux ratio depends on the type of source dominating 
the energetics of the interstellar medium. A mid-infrared diagnostic diagram first 
put forth by Genzel et al. (1998), and later explored by Peeters, Spoon, & Tie- 
lens (2004), plots the emission line ratio [O IV]25.89 /im/[Ne II]12.81 /im versus 
the strength of a mid-infrared PAH feature. In such plots AGN sources show en- 
hanced [O IV]25.89 fim emission and comparatively weak PAH feature strength. 
However, those two studies focussed on AGN-dominated and starburst nuclei. 
The upper panel of Figure ^ uses the 6.2 /jm PAH feature and the emission 
line ratio [O IV]25.89 jum/[Ne II]12.81 /im in a similar mid-infrared diagnostic 
diagram, but one that incorporates "normal" (starbursting/star-forming) nuclei 
and H II regions — the sensitivity of Spitzer allows us to probe to far fainter levels 
than heretofore possible. As expected, normal nuclei and H II regions extend 
the previously-observed trend: lower-luminosity star-forming nuclei and H II re- 
gions exhibit comparatively large 6.2 /im equivalent widths and relatively low 
ratios of [O IV]25.89 /im/[Ne II] 12.81 //m, indicating strong contributions from 
[Ne II] 12.81 /im cooling of H II regions and their PAH-rich photodissociation 
region surroundings, and negligible emission from AGN. 

High ionization lines like [O IV] 25.89 /im are somewhat difficult to detect 
in many SINGS sources. An alternative diagnostic diagram is provided in the 
lower panel of Figure ^ involving the more easily detectable [Si II] 34.82 /mi 
line, which has an ionization potential of 8.15 eV. As pointed out by Maloney, 
Hollenbach, & Tielens (1996), the [Si II]34.82 /im line is a strong coolant of X-ray 
irradiated gas. In X-ray dominated regions the X-ray emission dominates that 
from the comparatively small H II-like regions surrounding the hard-spectrum 
source. Moreover, X-ray dominated regions can be quite large since hard X-ray 
photons penetrate large column densities, and the conversion of X-ray energy to 
infrared continuum and line emission can be very efficient. Maloney, Hollenbach, 
& Tielens (1996) predict [Si II]34.82 /im to be prominent cooling line within X- 
ray dominated regions. We take advantage of this concept to present a new 
technique for distinguishing between AGN sources and star-forming regions. 
This technique relies on an easily-detectable, prominent cooling line of a low- 
ionization species associated with X-ray dominated regions, the dense interstellar 
material illuminated by power-law radiation fields. 

The dotted lines in Figure ^represent a variable mix of an AGN nucleus and 
a "pure" star-forming region. The 100% AGN anchor for these mixing models 
have values of EW(PAH 6.2 /im)= 0.01 /im, [O IV]25.89 /im/[Ne II] 12.81 /im=0.3, 
and [Si II]34.82 /tm/[Ne II] 12.81 /mi=2, whereas the 100% star formation an- 
chor has EW(PAH 6.2 /mi)~0.7 /im, [O IV]25.89 /iin/[Ne II]12.81 /iin=0.01, 
and [Si II]34.82 /im/[Ne II]12.81 /im=0.4. The dashed line in the upper panel 
of Figure ^ shows the approximate mixing model of Genzel et al. (1998; their 
Figure 5), obtained after empirically deriving a relation between their 7.7 /im 
PAH 'strength' (line-to-continuum ratio) and the 6.2 /im PAH equivalent width. 
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Many of the high ionization line data presented by Genzel et al. (1998) were 
upper limits, so it is unsurprising that their original curve lies above our curve. 
Short solid lines roughly perpendicular to the dotted AGN/star-forming curves 
delineate three regions in both panels of Figure ^ The boundaries are: 

Region[ I - II] : log([0 IV]25.89/mi/[Ne II]12.81/mi) = 10 log(EW[6.2/miPAH])+8.0 

(1) 
Regional - III] : log([0 IV]25.89/im/[Ne II]12.81/mi) = 3.41og(EW[6.2/zmPAH])-0.2 

(2) 
Region[IV - V] : log([Si II]34.82/mi/[Ne II]12.81/mi) = 5.01og(EW[6.2/miPAH])+4.8 

(3) 
Region[V - VI] : log([Si II]34.82/xm/[Ne II]12.81/mi) = 1.71og(EW[6.2/miPAH])+0.5 

(4) 
The population statistics for these regions show that Regions [I+IV] and [III+VI] 
are respectively representative (at the ~90% level) of Seyferts and star-forming 
systems. Regions [11+ V], on the other hand, contain a mix of classifications and 
thus represent transition regions — either the source classifications in this region 
are ambiguous or the region simply contains a more heterogeneous mixture of 
pure types. 

A Neon, Sulfur, and Silicon Diagnostic If the neon excitation is plotted as a 
function of [S III]33.48 /xm/[Si II]34.82 /an (Figure^), a more obvious separation 
of the star-forming and AGN-powered data points is observed. Not only do the 
low-metallicity Magellanic Cloud regions exhibit a higher neon excitation, nearly 
all of the "pure star-forming" nuclei and extra-nuclear regions show relatively 
elevated ratios in [S III]33.48 /mi/ [Si II]34.82 /mi. Note in addition that many 
of the filled squares representing starbursting/star-forming nuclei are located 
between the H II regions and the AGN. The boundaries are defined by curves 
with the same slope but differing offsets: 

log([Ne III]15.56/mi/[Ne II]12.81/mi) = 8.41og([S III]33.48/tm/[Si II]34.82/mi)+7, 

(5) 
where 7 =[+3. 3, +1.2, —2.5] for the lines demarcating Regions [I— II, II— 111,111— IV]. 

These results can be partially understood in the context of the cooling line 
physics introduced above. The [Si II] 34.82 /iin line is a significant coolant of 
X-ray ionized regions or photodissociation regions (Hollenbach & Tielens 1999), 
whereas the [S III] 33.48 /im line is a strong marker of H II regions. In other words, 
extra-nuclear regions and star-forming nuclei will show strong signatures of the 
Stromgren sphere coolant [S III]33.48 /mi, while AGN and their associated X- 
ray dominated regions or dense photodissociation regions will exhibit relatively 
strong [Si II]34.82 /tm emission in analogy to the increased strength of [O I] 6300A 
emission in AGN (e.g., Veilleux & Osterbrock 1987). In addition, the fraction 
of photodissociation regions falling within each beam will play a role in the 
line ratios. The data for Magellanic Cloud and Galactic H II regions stem 
from smaller physical apertures and thus are likely to have fractionally higher 
contributions from Stromgren spheres than photodissociation regions. 

Metallicity may be a factor as well. Since the central regions of galaxies 
typically are more abundant in heavy metals (Pagel & Edmunds 1981; Mc- 
Call 1982; Vila-Costas & Edmunds 1992; Pilyugin & Ferrini 1998; Henry k 
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Worthey 1999), and as explained above a lower metallicity can lead to harder 
radiation fields and thus enhanced high-ionization-to-low-ionization line ratios, 
it is possible that this AGN^H II nucleus^H II region sequencing along the 
[S III] 33.48 |Um/[Si II] 34.82 /im axis is affected by metallicity. However, the 
lower metallicity Magellanic Cloud data are not substantially to the right of 
the Galactic H II region data, so the effect is not solely due to metallicity. Al- 
ternatively, perhaps some of the star-forming nuclei have contributions from 
undetected weak AGN and thus are not "pure" star-forming nuclei, resulting 
in a location for star-forming nuclei on this diagram between AGN and H II 
regions. 



3. Summary 

We have presented mid-infrared diagnostic diagrams for a large portion of the 
SINGS sample supplemented by archival ISO and Spitzer data. A portion of our 
work solidifies and extends previous ISO-based mid-infrared work to lower lumi- 
nosity normal galaxy nuclei and H II regions using the Spitzer Space Telescope. 
We also present new diagnostics that effectively constrain a galaxy's dominant 
power source. The power of the diagnostic diagrams of Genzel et al. (1998; see 
also Peeters, Spoon, & Tielens 2004) for distinguishing between AGN and star- 
forming sources in dusty ULIRGs is that mid-infrared lines and PAH features 
are much less affected by extinction than their optical counterparts in a tradi- 
tional diagnostic diagram. Unlike diagrams put forth by Genzel et al., which 
rely on detecting relatively weak high- ionization lines like [O IV] 25.89 /im and 
[Ne V] 14.32 /im, we provide a new diagnostic that utilizes a strong /owMonization 
line. The advantage of using a line ratio like [Si II] /[Ne II] is that singly-ionized 
silicon and neon respectively have ionization potentials of only 8.15 and 12.8 eV, 
so they can both be observed over a large range of physical conditions. This is 
similar in concept to previous efforts that have taken advantage of [O I] lines (e.g., 
at 6300 A or 63 /im) that are coolants of the X-ray dominated regions (or dense 
photodissociation regions) surrounding AGN. In plots of both [O IV] /[Ne II] vs 
PAH equivalent width and [Si II] /[Ne II] vs. PAH equivalent width we identify 
a region where >90% of the sources are Seyfert or LINER. Likewise, another 
region in both plots is >90% H II regions or star-forming nuclei. 

Another useful mid-infrared diagnostic is [Ne III] 15.56 /im/[Ne II] 12.81 /im 
vs. [S III] 33.48 /im/[Si II] 34.82 /im. This plot tracks the excitation power of 
the radiation field on one axis, while the other axis is a relative measure of the 
cooling of H II regions and X-ray dominated regions (or dense photodissociation 
regions). Similar to what is found for the [Si II] /[Ne II] vs. PAH diagnos- 
tic alluded to above, both starbursting nuclei and extranuclear regions stand 
apart from nuclei that are powered by accretion-powered disks. Moreover, com- 
pared to starbursting nuclei, extranuclear regions typically separate even further 
from Seyfert nuclei, especially for low-metallicity environments. Presumably this 
extranuclear < — muclear separation occurs since extranuclear regions are cleaner 
representatives of H II regions than starburst nuclei, because their stellar popu- 
lations and interstellar medium structure are less complex. Extranuclear regions 
more likely contain younger stellar populations since they trace a single burst, 
as opposed to the average of multiple star formation episodes for nuclei (e.g., 
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Dale et al. 2004b). Finally, we note that it is difficult to clearly distinguish 
between pure Seyfert and pure LINER sources using these diagnostics. 
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Figure 1. The ratios of mid- infrared forbidden lines as a function of the 
6.2 /im PAH feature equivalent width. SINGS data are displayed with la- 
error bars based on the statistical uncertainties; archival data are displayed 
without error bars. The dotted lines are linear mixing models of a "pure" 
AGN and a "pure" star- forming source (see text). The dashed line in the top 
panel is a similar mixing model first presented by Genzel et al. (1998). The 
solid lines delineate regions distinguished by Seyferts, LINERs, star formation, 
etc. (see S l2~TT) . 
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Figure 2. A neon, sulfur, and silicon diagnostic diagram involving ratios of 
lines at different ionizations is displayed. The data are displayed as described 
in Figure n The lines delineate regions distinguished by Seyferts, LINERs, 
star formation, etc. (see § 12. 1.(1 . 



